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A new method is presented for Fourier decomposition of the Helmholtz 
Green Function in cylindrical coordinates, which is equivalent to obtain- 
ing the solution of the Helmholtz equation for a general ring source. The 
Fourier coefficients of the Helmholtz Green function are split into their half 
advanced+half retarded and half advanced— half retarded components. 
Closed form solutions are given for these components in terms of a Horn 
function and a Kampe de Feriet function, respectively. The systems of 
partial differential equations associated with these two-dimensional hy- 
pergeometric functions are used to construct a fourth-order ordinary dif- 
ferential equation which both components satisfy. A second fourth-order 
ordinary differential equation for the general Fourier coefncent is derived 
from an integral representation of the coefficient, and both differential 
equations are shown to be equivalent. Series solutions for the various 
Fourier coefficients are also given, mostly in terms of Legendre functions 
and Bessel/Hankel functions. These are derived from the closed form 
hypergeometric solutions or an integral representation, or both. Numer- 
ical calculations comparing different methods of calculating the Fourier 
coefficients are presented. 
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Abstract 



1 Introduction and overview 



The inhomogeneous Helmholtz wave equation is 




(1) 
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and this has the well known free-space retarded Green function [1, p. 284] 



G H 03,r-r') 



exp (iff | r — r' |) 



(2) 



47T | r — r' | 



where r is a field point, r' is a source point and (3 is the wave number, here 
considered to be a general complex number. The free-space Green function {2J 
is restricted to values of (5 such that | Gh (/3, r — r') |— > as | r — r' |— > oo. For 
general dispersive waves with (3 = a + ia where a and a are real, then a ^ 
is a condition for this to hold. In the limit as (3 — ► these equations reduce 
to the Poisson equation and its corresponding Green function Gp (r — r'). The 
general retarded solution <E> (/3, r) of the Hclmholtz equation at a field point 
r for a general source density p(r), subject to the boundary condition that 
$ (/?, r) — ► as | r |— > oo, is given in terms of the Green function as 



where the volume integral is to be taken over all regions of space where the 
source density p (r) is non zero. 

Many problems of practical interest have some element of axial symmetry 
and are best treated in cylindrical coordinates (r, </), z), the Cartesian compo- 
nents (x, y, z) of r being related to the cylindrical components by (x,y,z) — 
(r cos <fi, r sine/), z). It follows immediately from this relation that the distance 
between a source point and a field point is given by 



The solution for $ (/3,r) when p(r) is a general circular ring source is of par- 
ticular interest, with applications such as circular loop antennas [2], [3], [4], the 
acoustics of rotating machinery [5] and acoustic and electromagnetic scattering 
[6]. For the simpler Poisson equation most of the analytical solutions found in 
the literature for cylindrical geometry are either ring source solutions or can 
be easily constructed from them by integration or summation. Examples are 
gravitating rings and disks, ring vortices and vortex disks, and circular current 
loops and solenoids. 

The source density p c (r, R, z) for a thin circular ring of radius R located in 
the plane z — Z is of the form 



where / (<f)) is the angular distribution of the source strength around the ring. 
This can be most conveniently described by a Fourier series of the form 




(3) 



| r — r' |= J r 2 + r' 2 + (z — z') 2 — 2rr' cos (0 — <fi') 



(4) 



Pc (r,R, Z) = f{4>)8{r-R)5{z-Z) 



(5) 



OO 



a m cos (rrufi) + b m sin (mcf))) . 



(6) 
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where the Fourier coefficients a m and b m are given by [7, p. 1066] 

a,n = 2tt I f ^ C ° S ^ m< ^ d ^ ^ 



o 

2;r 



b m = ^- //(</») sin (m£)#. (N) 

Z7T 
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From equation (j4]), the Green function (|3|) is even in the variable ip = <j> — <j>, 
where <f> is the angular coordinate of r and <j> is the angular coordinate of r'. It 
is convenient to exploit this symmetry when substituting equations ((5]) and ([6]) 
into From the identity / (<fi + ip) = / (0') we obtain 



/ (^o = -^ + ^2 (°™ cos ( m< ^) + b ™ sin ( m< / ) )) c ° s 

m— 1 



+ (— a m sin (m<j>) + b m cos (mcj))) sin (mi/i) (9) 

and on substituting equations ^ and ^ into an d performing the volume 
integration, the odd terms proportional to sin (mip) in equation ([9]) do not con- 
tribute to the solution $ (/3,r) as Gh (/3, r — r') is even in The remaining 
integrals from the even terms can be calculated over the reduced interval from 
to 7r. This gives the solution <£> c (j3, r,R, Z) of the Helmholtz equation for a 
circular ring source with general / (0) in the form 

<f c (/3, r, R,Z)=-^G° H (/3,r,R,z- Z) 

OO 

- ^2 cos ("V) + b ™ sin (™4>)) g h (0) r,R,z- Z) (10) 

where 



* „exp I i(3\J r 2 + R 2 + (z - Z) - 2ri?cos ip 

G% (ft, r,R,z-Z) = - / ^ — L cos (mip) dif) 

n { yV 2 + R 2 + (z - Zf - 2rRcost/j 

(11) 

and where the explicit dependence of the solution on the constant ring parame- 
ters R and Z has been introduced in these definitions. Introducing the Neumann 
factor e m such that e m = 1 for to = and e m = 2 for m > 0, and defining bo = 
allows ([TUl) to be expressed more concisely as 



2 

m— 



- oo 

$ c (/?, r, i?, Z) = -- Om cos (mcj)) + b m sin (to0)) e m G% (/?, r,R,z-Z). 

(12) 
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Apart from a constant factor, the terms e m G^ (/?, r,R,z— Z) in (112[) are also 
the coefficients in the Fourier expansion of the Green function ([2]) itself, when 
the source point is given by r' = (i?, Z). From equations ©, and © 
this is given by 

oo 

G H (P,r-r') = -—J2 ^Gh (ft r,R,z- Z) cos (m (0 -<£')). (13) 

m— 

Thus the solution $ c (/3, r, 72, Z) of the Helmholtz equation for a general ring 
source can be constructed directly from the coefficients G]} (/?, r,R, Z — z) in 
the Fourier expansion of the Green function ([3]) . This provides in large measure 
the motivation to analytically construct the Fourier series for the Helmholtz 
Green function. 

For the Poisson equation with (3 = the corresponding Fourier expansion of 
the Green function has already been given in closed form as [8] : 

1 °° 

G P (r - r') = - V e m Q m „ 1/2 (w) cos (m (0 - 0')) (14) 

47rVrJZ ^ 

where 

.= - 2 + J?2 + ^-^ (15) 

is a toroidal variable such that u> ^? 1 and the Q m ~i/2 (<*>) ar e the Legendre 
functions of the second kind and half-integral degree, which are also toroidal 
harmonics. The Fourier expansion given by equations (| 14[) and (|15[) can be 
obtained immediately by writing the Green function (|2|) for (3 = in the form 

G P (r - r') = 1 = (16) 



&irVrRy/2w - 2 cos (0 - <j>') 

where u> is given by (115|) . and noting that the function Q rn -\/2 (<- L has the simple 
integral representation [7, eqn 8.713] 

/ cos (mjj) di/j 
7 v 2w — 2 cos ^ 



An alternative derivation of (fl"4"l) employs the Lipschitz integral [7, eqn 6.611 1] 

oo 

f J {sa)exp(-s\b\)ds=-=±= (18) 
o 

and Neumann's addition theorem [9, eqnll.2 1] 

oo 

Jo [s \Jr 2 + R? — 2Rr cos ipj = e m cos (m"0) J rn (sr) J m (sR) (19) 



m=0 
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to obtain the well known eigenfunction expansion 

G P (r-r') = -±-x 

oo 00 

e m cos (m (0 - </>')) / J m (sr) J m (sR) exp (-s | z - Z |) ds. (20) 

m=0 q 

This reduces to (fT4"f on employing the integral [7, eqn 6.612 3], [9, eqn 13.22]: 

00 

J m (sr) J m (sR) exp (-a | Z - z |) rfs = — i=Q m -i/2 M ■ ( 21 ) 



The generalization of (|20|) for the Helmholtz case is also well known [10, p. 888] 

Gn(Ar-r') = ~x 



00 00 

^ e m cos (m (0 - 0')) / exp (i\Z — z\ \J 1 - s 2 ) J m (sr) J m (sR) 

m=0 VP 

(22) 

This can be similarly obtained from Neumann's theorem by employing the in- 
tegral [7, eqn 6.616 2] 

exp (-ax) J Q (bVx* - l) dx = ex P (~^ + 6 1 ( 23 ) 
1 

instead of the Lipschitz integral. Equation (f2"2")) gives the Fourier coefficients of 
the Helmholtz Green function in the form 

00 

G£ ((3,r,R,z-Z) = ijexp(i\Z-z\ \j p 2 - s^j J m ( sr ) J m (sR) 


This reduces to (|2"0|) in the limit as (3 — * but unfortunately the integral in 
is not given in standard tables for (3 7^ 0. Numerical evaluation of this integral 
requires care, as the integrand is oscillatory and singular in an infinite range 
of integration, though the integrand tends exponentially to zero as s — > 00. 
Equation (jlip is a convenient alternative numerical evaluation of the Fourier 
coefficients, provided m is not too large. 

The integrals (TTTj) and (f2"4"| contain the additional parameter f3 which is not 
contained in (fTTj) and (|2"Tj) . As a consequence of this, the closed form general- 
ization of (|14|) for the Helmholtz case involves two-multidimensional Gaussian 
hypergeometric series, and the main purpose of this article is to present these 



(24) 
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solutions and various related results. The core idea leading to the solution is 
expansion of the exponential in ([2]) as the absolutely convergent power series [4] 



n=0 

where 

| r-r' | n_1 = (r 2 +R 2 + (z- Z) 2 -2rR cos 4>) . (26) 



Hence 

G£ (0, r,R,z-Z) = y2 ^-Im.n (r, R,z~Z) (27) 

n=0 

where 

7T 

If/ o \ (™-l)/2 

I m ,n (r, R, z - Z) = - / (r 2 + i? 2 + (z-Z) -2ri?cos'0) cos (mip) dtp . 

o 

(28) 

The integral (|28[) can be evaluated as a series by binomial expansion and this 
gives a double series for the Fourier coefficient G^ (0, r, R, z — Z). The expan- 
sion of (|28[) gives an infinite number of terms for n even and a finite number of 
terms for n odd. These two cases are best treated separately and it is therefore 
convenient to split the summation over n in (|25p into odd and even terms. This 
is equivalent to splitting the Green function ^ such that 

G H (0, r - r') = A+ (0, r - r') + A_ (0, r - r') (29) 

where 

A (R r. ,'\ 1 f exp(i/?|r-r'|) , cxp (-»/3 | r - r' |) \ 

A + (Ar- r ) = -_ ^ _ — + _ — j (30) 

is the half advanced+half retarded Green function and 

A (a„ j\ 1 ( exp(z/3 I r- r' |) exp(-z/3| r-r' |) ^ 

A - (/3 < r - r } = - ^ (, — — J (31) 

is the half advanced— half retarded Green function. The corresponding Fourier 
coefficients are split in the same manner such that 

A™ (0, r,R,z-Z) = ± (Gh (0, r,R,z-Z) + G% (-0, r,R,z- Z)) (32) 

K m (0, r,R,z-Z) = ± (Gfj (0, r,R,z-Z) + G£ (-0, r,R,z- Z)) . (33) 

For real 0, splitting the Green function in this way is equivalent to dividing it 
into its real and imaginary parts, but this is not the case for general complex 0. 
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It is shown in Section 2 that the Fourier coefficients in (|32p and ([55)) are given 
respectively by 

A™ (rR : 7, k) = 

2m+l 



and 



T(m + 1/2) (k 
mW itrR \2 

A™(ri?, 7 ,fc) 



H 3 m+ 1/2, m+ 1/2,2™ + l,k 



2 7 



(34) 



jk 



2m+l 



0:1:0 
1:1:0 



VF(2m+ 1)! V 2 
m + 1/2; -; j 2 k 2 



m + 3/2: 2m + 1; 



where 



4ri? 



(r + i?) 2 + (z - Zf 



and 



7 = /?V^ + ^) 2 + (^-^) 2 



A± (ri?, 7 ,fc) = A£ (/3,r,i?,z-Z) 



(35) 

(36) 

(37) 
(38) 



The variable k is the usual modulus contained in elliptic integral solutions of 
elementary ring problems and is related to the toroidal variable lo by 



2-fc 2 
k 2 



(39) 



The function H3 in equation (|34[) is one of the standard Horn functions [11, eqn 
5.7.1 31] and is equivalent to the double hypergeometric series 



H 3 ( m + 1/2, m + 1/2, 2m + 1, k 2 , 1- 



EE 

n=0 p=0 



(m+l/2)„_ p (m+l/2) n , ^_ 

(2m + l)„n!p! ^ ' ' " 



2 \ p 



(40) 



The Kampe de Feriet function [12, p. 27] in ([33)) is equivalent to the double 
hypergeometric series 



F 0:l;0 
r 1:1:0 



- : m + 1/2; -; j 2 k 2 -7 2 " 1 
m + 3/2: 2m+l; - 



(m + l/2) n 



SS( TO + 3 / 2 Wf (2m + l)„n!p! 



4 ' 4 

2 ; 2 \ n 

Yk 



(41) 
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The integral (|2"8")l can also be evaluated using an integral representation for 
the associated Legendrc function of the first kind, and it is shown in Appendix 
A that this gives the series expansion: 



(w 2 - l) 1/4 V2ri? 



,2 



1 



l/4\" 



r vj >_L r((n + i)/2) , u, , 

^ n! r(m+(n + l)/2) (n_1)/2 Iv/^ ' 

n— x 



where 



A = /3V27R (43) 

and 

G^{rR,\u) = G^{p,r,R,z- Z) . (44) 

The Legendre function in equation (|4*^|) reduces to an associated Legendre poly- 
nomial for odd n. The series in (|4"2"]l can be split into even and odd terms such 
that 

Gh {rR, A, u) = ^ (A™ (ri?, A, fc) + A™ (riZ, A, w)) (45) 

where 

A^(r J R,A, W ) = A^(/3,r,i?,z-Z) (46) 

and it is shown in Appendix A that the even and odd series can be expressed 
respectively as: 

4 m / „ , (-lT^/ AVg^i y <&-i/ 2 M 

+ lr ' ,LJ) yfrR tl\ 4 j p!r(p-m+l/2)r(p + m + l/2) 

(47) 



p=0 



A-(r J R,A,o;) = i 7 4-^ 



(-ir ^ ( \ 2 ^—i\ p+m+1/2 qT.T h 



~^ V 4 / p!r(p + m + 3/2)r(p+ 2m + 1)' 

(48) 

In equation (|48[) the Legendre function is purely imaginary for real A. In the 
static limit as A — > then A™ (ri?, 0, w) = and from the gamma function 
identity [7, eqn 8.334 2] 

r(l/2-m)r(l/2 + m) = (-l) ro 7r [m G No] (49) 

then equation (|4"7| reduces to 

Ay(riz,o,M) = [ meNo ] (50) 

ny/rR 



as it must do for consistency with ([14]) . 

The solutions in terms of two-dimensional hypergeometric functions defined 
by Equations (|34 |) -([38 |) and (|40|) -(|4T |) can be summed over either index to give 
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the solutions as series of special functions. It is shown in Section 3 that summa- 
tion over the index n in equation ([40]) gives equation (|47|) . exactly as given by 
the integral representation. However, summation over the index n in equation 
([IT]) gives instead the series solution 



A™ (T\R,7,fc) 



m\\frR \ 4 



, \ 2m+l 
7K 



e r(P+m+ 3/2)p! ) ^ r + 1/2; n+m+ 3A 2m + 1: V 

(51) 

A hypergeometric identity to reduce the hypergeometric function in equation 
(|5ip to other well-known special functions does not seem to be available in stan- 
dard tabulations. It might nevertheless be conjectured that l[5"T]) could somehow 
be reducible to equation ([4"8|) , but this is not in fact the case. It is easily ver- 
ified numerically that although equations (|4"5)) and ([51]) both converge rapidly 
to the same limit, the individual terms do not match. Hence, equation (|51[) is a 
distinct series from equation (I48[) . It is also shown in Section 3 that summation 
over the index p in equations (I40[) and (|4"Tj) gives the Bessel function series: 

Xmr u i \ -1 f r(n + m+l/2) / 7 P\" +m+1/2 v 

A? (rR, 7, k) = — = > — -t-^ - — Y„, m , i (7) 



n=Q 

(52) 

A-(,n 7 ,*) ■ g-r(,. + n , + 1 /2) /^y— » 

~ l j 2\A 7 R Z - I r n + 2m+l n! V 2 / n+m+ 3 w; 

v n— / \ ✓ 

(53) 

and these two series can be conveniently combined to give a series of Hankcl 
functions of the first kind: 

Gg (rfl, 7, fc) = -!= Y T n + m + ffi f ^ ff« . (7 ). 

2\/r r R ^ T (n + 2m + 1) n\ V 2 / n +™+i w/ 

v n— / \ / 

(54) 

From the solutions ([34")) and ([33)) it can be seen that dimensionless Fourier 
coefficients defined by 

5 £(fc 2 , 7 2 /4) = VTRK^ (ri?, 7 ,fc) (55) 

depend only on the two dimensionless variables x = k 2 and j/ = 7 2 /4. The 
functions <?™ (x,y) are given explicitly by equations ([34]) and l|35p as: 

9^ (x, y) = y£±±^ x a H a (a, a, 2a, x, y) (56) 



and 



q + 1: 2a: -; Xy '~ y 



(57) 
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where 

a = m+l/2. (58) 

Two-dimensional hypergeometric series such as (|56[) and (|57[) are associated with 
pairs of partial differential equations [11, section 5.9] and these can be used to 
construct ordinary differential equations for g± (x, y) with y fixed and x as the 
independent variable. It is shown in Section 4 that for constant y the coefficients 
g± (x,y) both satisfy the same fourth-order ordinary differential equation in x: 



(1 — x) x 



dx A 



(6 - 9x) x 



: a 9± 
dx 2 



dx s 



(a(l 



6 — I8x — xy + 2y) x 



dgl 



(2a (1 - a) - 2x (3 - y)) x^-+(y 2 + 2a (1 - a) y - 3a (a + 1) (a + 2)) g% = 0. 

(59) 



dx 

In Section 5 an integral representation is derived for 
y m (\,Lo) = irV2^RG% (rR,X,cj) 



(60) 



and this is used to derive a fourth-order ordinary differential equation for y m (A, to) 
in terms of lo: 



(I-- 2 ) 



diu 4 



r . d y m . 2 
-6a; , „ +(m z - 
duj 6 



X 2 u 25 d 2 y m 
2 4 ' diJ 1 



>. dy m 
du> 



y m = o. (61) 



In the static limit as A — > 0, equation (|6T|) reduces to 



dw 2 



1-^ 



d 2 y„ 
dw 2 



- 2uj 



dy m 
dw 



m 

4 



= 



where 



d 2 y n 
duj 2 



2to 



dy m 
dto 



Un 







(62) 



(63) 



is Legendre's equation of degree m— 1/2 [7, eqn 8.820]. It is also shown in Section 
5 that the differential equations (|59| and (folj) . obtained by quite different routes, 
are equivalent. The special functions used in the analysis are given in Table 1. 

Recurrence relations for the Fourier coefficients for the Hclmholtz equation 
were investigated by Matviyenko [6], but the closed form solutions and differ- 
ential equations presented here appear to be new. Werner [3] presented an 
expansion of the Fourier coefficient as a series of spherical Hankcl functions, su- 
perficially similar to equation (|54p . but the two expansions are distinct. The two- 
dimensional hypergeometric series approach applied here to obtain the Fourier 
expansion for the Helmholtz Green function has recently been applied to obtain 
the Fourier expansion in terms of the amplitude 4> for the Legendre incomplete 
elliptic integral of the third kind [13]. 

The numerical performance of the various expressions for the Fourier co- 

® 



efficients was investigated using Mathematica 
Appendix C. 



[14] and this is examined in 
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F 0:l;0 
r l:l:0 



Symbol 

(«)„ 
B(x,y) 

2-Fi (a, 6; c; x) 
1^2(0; b, c; a;) 
- : a; -; 
6: c; -; 
HP (x) 
H 3 (a, b,c,x,y) 
J v (x) 
P?(x) 

Y v {x) 
Y{x) 
8{x) 



Special Function 
Pochammer symbol 
Beta function 
Gauss hypergeometric function 
A hypergeometric function 

x, y A Kampe de Feriet function 

Hankel function of the first kind 
The H3 confluent Horn function 
Bessel function of the first kind 
Associated Legendre function of the first kind 
Associated Legendre function of the second kind 
Bessel function of the second kind 
Gamma function 
Dirac delta function 

Table 1: Special Functions Used 



2 Solution in terms of two-dimensional hyper- 
geometric series 

The power series expansion (|2~T1) for the Fourier coefficient G^ (/3, r,R, z — Z) 
can be expressed in the form 



es(r*, 7 ,*) = t^£M^ 



7T/2 



7rv rR 



n=0 



nl 



cos(2m#)(l-fc 2 sin 2 #) ( "- 1)/2 



d9 (64) 



where 



G% ((3, r,R,z-Z) = G m (rR, 7 (p, r,R,z-Z),k (r, R,z- Z)) 



(65) 



and k and 7 are defined by ([55 ]) and ([57 ]) . The term (l - /c 2 sin 2 0) ( " 1)/2 in (M]) 
can be expanded binomially to give G^ (ri?, 7, k) as a double series containing 
integrals of the form 



tt/2 

T m!P = J sin 2p cos (2m6>) d<9. 



(66) 



This integral is given by Gradshteyn and Ryzhik [7, eqns 3.631 8,12] in a form 
which can be recast as 



j (-1)'"* 1 

m ' p 2 2 P+ 1 (2p+l)B(p + m+l,p-m+l) 

Imp — for p < m 



for p m (67) 
(68) 
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and expressing the beta function in (|67p in terms of gamma functions and em- 
ploying the duplication theorem [7, eqn 8.335 1] 

2 2x_1 

r (2x) = —^r (x) r ( x + 1/2) (69) 

V 71 " 

gives after some reduction the alternative form 

(-i) m yi r(p + i/2) P ! 

l m v = ; ; ; r for p >■ m. 70 

' p 2 r(p + m+l)rO-m + l) F V 7 

The binomial expansion of (|64p gives an infinite series for n zero or even and 
a finite sum for n odd, and these two cases must be treated separately. It is 
therefore convenient to split the series for such that 

G£ (ri?, 7 , k) = A™ (ri?, 7, fc) + A m (ri?, 7 , fc) (71) 

and on employing equation l|69p the divided series are given by 

(-!)«*- (- 7 2 /4) n <wn^ ,2^2M«-^ 



A™ (ri?, 7 , fc) = V Y ) ' ' ,\ , / cos (2m9) (l - fc 2 sin 2 0)"~ 1/2 dfl 



L 7n,s 



,2 (?2) 

A m (ri?, 7, fc) = - ( ~ X )^: 7 y (~ 72/4 \ , / cos (2m9) (l - fc 2 sin 2 9)" d9. 
- K ' 2V^R „t^r(n + 3/2)n! 7 V j 

(73) 

Binomial expansion of the integrals in equations (|72| and (f73|) gives respectively 

(- 72/4 )" r(^n + i/ 2) (fc 2 ) s 

A +l ri ™*>- ^fl^rln + W r(l/2-n) .! 

(74) 

L . (r47lk)a ffl^M^ ( 75 ) 
2\fwrR n r (n + 3/2) si (n - s)\ 

and employing the explicit formula (|70[) for i m a in (jT4"|) and ([75]) gives respec- 
tively 

A"(ri?, 7 ,fc) = — 

» » r(.s-n+l/2)r( S + l/2) ( 7 2 /4)"(fc 2 ) a 
^ ^ r(s + m + l)r(s-m + l) n! 1 ' 

n— s—m x / ' 

A™ ( rjR)7 ,A ; ) = -^ 7 =x 
4Vri? 

f> A r(.s + i/2)(- 7 2 /4)"(-fcr 

4^ 4^ ( n - s ) !r + 3 /2) r (s + m + 1) r (s - m + 1) 1 j 



n— m s—m 



where the gamma identity (|49p has been used to simplify equation (|76p . 



12 



where in ([78 



2.1 The series for A™ (rR, 7, k) 

The substitution s = p + m in equation (|76p yields after some reduction the 
double series: 

mWirrU \ 
- » (m + l/2) p _ n (m+l/2) p (fc 2 ) P ( 7 2 /4)" 

(a), S ^ (79) 

is the Pochhammer symbol. The double hypergeometric function in (|78p can be 
identified as one of the confluent Horn functions [11, eqn 5.7.1 31] and hence 

Am, D M r(m+l/2) /fc\ 2 " 1+1 

H 3 ^m+l/2,m+l/2,2m+l,/c 2 ,y^ . (80) 

The convergence condition given in [11, eqn 5.7.1 31] for the double series in 
equation (|T8[) is fc 2 < 1, which always holds. The order of summation in equation 
l|78p can be reversed, but the order of the arguments k 2 and 7 2 /4 in equation 
([50)) cannot be exchanged. 

2.2 The series for A m (rR, 7, k) 

Equation ([77)) can be converted to a doubly infinite series by reversing the order 
of summation, which gives 

A™ (rR n ,k) = -!%x 
WrR 



00 00 



y y r (s + i/2)H 2 /4)"(-fc 2 )- 

^ ^ (n-s)!r(n + 3/2)r(s + m + l)r(s-m + l)' ^ J 



s—m n—s 



The substitution n = s -\-p in (|5T|) gives 



A™ (r J R, 7 ,fc) = -^x 



00 00 



v v r (^ + i/2) (f#y (_f\ p 

£^p!r(p + s + 3/2)r( a + m + l)r( a -m + l) V 4 / \ 4 / 
and the further substitution s = ri + rn in Ij82|) gives 

^•^'VfT V) * 
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oo oc 

EE 



T(n + m+ 1/2) 



- ,A ^T(p + n + m + 3/2) r fn + 2m + 1) n!p! 

p— o v / / v / 



1 / 7 2 fc 2 



2\ P 



(83) 



Expressing equation (|83p in terms of Pochhammer symbols gives after some 
reduction the double hypergeometric series 



A™ (ri?, 7 ,fc) 



7 A: 



2m+l 



/t\R(2to+1)! V 2 
(m + l/2) n 1 fj 2 k 2 



2\ V 



(84) 



and this can be expressed as a Kampe de Feriet function as defined by Srivastava 
and Karlsson [12, p. 27]: 



pp:q;r 
l:m:n 



{a P } ■ {b q } ; {c r } ; 

W: {/U; hJ; ' y 



oo oo nw^nmiw 

u=l u=l u=l 



x 3 y l 



j\ i\ 



(85) 



i=o i=o n c«-)i-H n (/u- n (-a 

u— 1 u— 1 u— 1 

In the definition (|85p. {%,} = cji, ..a p and {b q } and so on, are the lists of the 
arguments of the Pochhammer symbols of the various types which appear in 
the products on the right-hand side of the equation. If a list has no members, 
it is represented by a hyphen. Comparing with (|55|) gives 



A m (rR,j,k) = 



F 0:l;0 
r l:l;0 



<rR(2m + 1)! 



1/2; -; 7 2 fc 2 - 7 2 



(86) 



m + 3/2: 2m + 1; -; 4 ' 4 

3 Consequences of the hypergeometric formulas 

In the static limit as (3 — ► then equation ([34]) reduces to 



A" (riZ, 0, jfe) 



r(m+ 1/2) /fc 



2m+l 



and from (|37p and the standard hypergeometric identity [15, eqn 7.3.1 71] 



2 Fi (m + 1/2, m + 1/2; 2m + 1; fc 2 ) 

(87) 



2 F X (a, 6; 26; z 



(88) 
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this reduces to 

1 „ /2-fc 2 



A+ (rJZ, 0, fc) - _L=Q m _ 1/2 ( (89) 
in agreement with equations (fT4| and ((37)) . 



3.1 Fourier coefficients as series of special functions 

The double series given by equation (|34[) can be summed with respect to either 
the index n or the index p in the definition (|40ll . Summing with respect to n in 
(|40|) gives a series of Bessel functions of the second kind: 



(90) 

where the gamma function identity (|49p and the Bessel function identity 

Y v (7) = [cos (j/tt) - J- v ( 7 )] (91) 

sin z/7T 



have been employed to obtain equation (|90|) . Summing instead over the index 
p in l|40p gives the alternative series 



A+ (ri?, 7 ,fc) 



1 /fc^ 2m+1 



7/7 



oo 12 I A\ n 

V r (to - n + 1/2) ^ 7 , ' 2 Fi (m+ 1/2 - n,m+ 1/2; 2to + 1; fc 2 ) . (92) 

n = 

This can be reduced using and (|4*9")l to give: 

~ /tVT^\ " Q"^i /2 ((2-fc 2 ) /fc 2 ) 

4 J n!r(n-TO+l/2)T(n + TO+l/2)' { ' 

The dimensionless variables 7 and A defined by equations ([37)1 and (|43|) respec- 
tively are related by 

7 = ^ (94) 



and substituting this equation and equation (|39j) in equation (|93|) gives imme- 
diately equation (|47)) . 

Summing with respect to p in equation (|83p gives the Bessel series 

A_ (ri?, 7, fc) - ^ r (p + 2m + 1) p! 2 J J ™ + H. ™ ^ 
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The corresponding summation over the index n gives 
A™ (rR, 7 ,k) = 



■ i— / i \ 2m+l 

l V 7T 1 1^ 



in 



WrR 



^ 1 f-l 2 \ n ( , , l 2 k 2 

E wn^^ffln. T iF 9 m+l/2;» + m + 3/2,2m+l- 



^ T (n + to + 3/2) n! V 4/1 '4 

n— f / \ / \ 

(96) 

There seems to be no hypergeometric transformation listed in standard tables 
suitable for directly reducing the hypergeometric function in this equation. 

Equations (|90|) and (|95|) can be conveniently combined to give a series of 
Hankel functions of the first kind: 

Qm (rRjk) -^f T( P + m+ l/2) (Jlf\ ^ 2 (1) 



where 



(i) 



(97) 



Kp (7) = J,n + ±+p (7) + ^m+i+p (7)- 



4 Differential equations from the two-dimensional 
hypergeometric solutions 

Erdelyi et. al. [11, section 5.9] tabulate the partial differential equations satisfied 
by all the functions in Horn's list. They employ the notation reproduced below 
for the various partial derivatives 

dz dz d 2 z d 2 z ^ d 2 z 

^ dx' dy dx 2 ' dxdy ' dy 2 

where z is any function on the list. Each function in the list satisfies two 
partial differential equations, and unfortunately those given in [11, 5.9 34] for 
z (a, /3,5,x,y) = H3 (a, j3, 8, x, y) contain typographical errors. The correct 
equations can be shown by the methods given in [11, section 5.7] to be: 

x (1 - x) r + xys + [5 - (a + [3 + 1) x] p + /3yq - a[3z = (98) 

yt-xs + (l-a)q + z = (99) 

and for the particular case considered here we have /3 = a and <5 = 2a so that 
(|98|) reduces to 

x (1 — x) r + xys + [2a — (2a + 1) x]p + ayq — a 2 z = 0. (100) 

Similar equations can also be derived for the Kampe de Feriet function defined 
by equation (l4Tj) . For the definition 



-I \ zti0:1;0 



- ■ <* - 
a + 1 : 2a: — ; 



(101) 
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and the notation 

dz dz d 2 z _ d 2 z - d 2 z 

OU OV OU z OUOV OV* 

then the equations corresponding to (|M)l and (| 1 00[) can be shown to be 

z = vi+us+ (a + l)q (102) 
z = lap + uf + q + vi. (103) 

4.1 Fourth-order differential equation from the Horn Func- 
tion H 3 (a, ft, 6, x, y) 

Writing 

z(a,x, y) = H 3 (a,a,2a,x,y) (104) 
then (|56]) can be expressed as 

o2ra+l | 

z (a, x, y) = - J x-g- (x, y) . (105) 
1 (m + 1/2) 

The ordinary differential equation for (x, y) in terms of a; = k 2 can be derived 
by first obtaining the corresponding differential equation for z (a, x, y) from (|99[) 
and poop and then substituting Q105P into this equation. Although straight- 
forward in principle, this procedure is rather intricate in practice, and only the 
essential elements of the derivation are given below. 

4.1.1 Differential equation for z(a, x, y) 

It is convenient to define a differential operator D such that: 

D = x4- (106) 

ax 

which has by definition the properties: 

Dx = x (107) 

Dy = -y (108) 

D(xy) = (109) 

Dz = xp — yq (HO) 

Dp = xr — ys (HI) 

Dq = xs-yt. (112) 
Combining equations (|99)) and (|112[) gives the equation 

Dq = z+(l-a)q (113) 
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and applying the operators d/dx and d/dy to this equation gives respectively: 

Ds=p-as (114) 

Dt = q + (2-a)t. (115) 
Eliminating the variable r gives a system of four coupled equations: 

Dz + yq = xp (H6) 

(1 — x) xDp = —xys + [(2a + 1) x — 2a] xp — axyq + a 2 xz (H7) 

Dq = z + (l-a)q (118) 

xDs — xp — axs. (H9) 
Eliminating p gives a system of 3 equations: 

D 2 z — xD 2 z+(2a — 1) Dz — 2oexDz+yz—xyz—a 2 xz = —xys+(l — a) yq (120) 

Dq = z + (l-a)q (121) 
xDs = Dz + yq — otxs. (122) 
Eliminating s gives the two equations: 

D 3 z - xD 3 z - (3a + 1) xD 2 z + (3a - 1) D 2 z -a(3a + 2) xDz + 2yDz 

+ a (2a — 1) Dz — xyDz — a 2 (a + 1) xz + 2 (a — 1) yz — axyz = —y 2 q (123) 

Dq = z + (l-a)q. (124) 
Eliminating q gives finally the fourth-order equation: 

(1 - x) D 4 z + [4a - (4a + 3) x\ D 3 z 

+ [5a 2 + a - 1 - (6a 2 + 9a + 2) x - xy + 2y] D 2 z 
+ [a (2a - 1) (a + 1) - a (4a 2 + 9a + 4) x + 2 (2a - 1) y - (2a + 1) xy] Dz 
+ [2a (a-l)y-a 2 (a + l)(a + 2)x-a(a + l)xy + y 2 ] z = 0. (125) 
This equation can be converted to standard differential form using the identity: 



which gives: 



Dz = x^- (126) 
ax 



d 4 z „ x ,. , ^d 3 z 



(1 ~ X) x4 dx^ + [2 {2a + 3) ~ (4a + 9) X] X "d^ 



+ [(a + 2) (5a + 3) - 3 (a + 2) (2a + 3) x - xy + 2y] x 2 — 
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dz 

+ [2a (a + 1) (a + 2) - (a + 1) (a + 2) (4a + 3) x + Aay - 2 (a + 1) xy] x— 

ax 

+ [2a (a - 1) y - a 2 (a + 1) (a + 2) x - a (a + 1) xy + y 2 ] z = Q. (127) 
From equation (|105p . the Fourier coefficient (x, y) is related to z (a, x, y) by 

z = Cx- a g™(x,y) (128) 

where the constant C is given by equation (| 1 05[) . Differentiating equation (|128p 
gives the relations 

dv dn m 

x ™=Cx- a [-ag™+x^} (129) 
ax ' T ax 



X 



= Cx~ a [a {a + 1) 5 ™ - 2ax^+ + x 2 ^+ ] 
dx^ efx dx^ 



3 d 3 z _ _ Q 

X ~~ : — r" — OX X 

dor 



do™ d 2 a m d 3 a m 

« (a + 1) (a + 2) 0™ + 3a (a + 1) x-|±- - 3ax 2 — ^- + ^~^\ (131) 



x 4 



: Cx- a X 



d A z 
dx 1 

dg™ 



[a (a + 1) (a + 2) (a + 3) g™ - 4a (a + 1) (a + 2) x-f :L + 

dx 

6a a + 1 x 2 — ±f-4ax 3 — ±f+x 4 -^ . 132 
dx z dx 6 ax 4 

and substituting these relations into equation (|36p gives after much reduction 
the fourth-order differential equation: 

1 - x x 4 — ^- + (6 - 9x x 3 — ?±- + 6 - a (a - 1) - 18x + y (2 - x x 2 — 

dx 4 dx J dx^ 



da m 

2 [a (a - 1) + x (y + 3)] x-gt + „a ff m = Q 



(133) 



4.2 Fourth-order differential equation from the Kampe de 
Feriet function 

Equation (|57|) can be written in the form 

9 m y) = l —rr. u a z (a, u, v) (134) 
(2m + ly. 

where 

u = xy (135) 

and 

v = -y. (136) 
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The equation for g™ [x, y) can be established in the same manner as for equation 
(|133|) , but having already derived (|133|) , it is enough to establish that g™ (x, y) 
and g™ (x, y) both obey the same differential equation. From equations (|135p 
and (fPgrJ)) then 

d_ _ ld_ 

du y dx 

d x d d 
dv y dx dy 



(137) 
(138) 



and with the definition 



w(a,x,y) =z(a,xy, -y) (139) 



and the notation 



dw _ dw _ d 2 w 

p = o - ; ? = q-; r = -5-9; s = ^r^r; * = 



9x ' <9y ' 9a; 2 ' dxdy ' <9y 2 
then we have 



p=- (140) 



2/ 

j 

y 



q=~p-q (141) 



f = ^ (142) 

2r 

5 = + 4 - - (143) 

y Z yi y 

2 

— ~ «Z/ iZv „. tZ/ _ - , , v 

t = t + 2—p+—r-2-s. (144) 

2r 2r y 

Employing equations (|139|) - (|144p in the Kampe de Feriet differential equations 
1(102)1 and |[T03]) gives 

yu) = — y 2 i + xys + crap — (a + 1) ycj (145) 

yw = (2a — x) p + x (1 — x) f — yc/ — y 2 i + 2xys. (146) 

To compare these equations with the partial differential equations from the Horn 
function, we note that 

g r ?(a,x,y) = Cx a z(a,x,y) (147) 

whereas 

g™(a,x,y) = Ex a y a w(a,x,y) (148) 
where C and E are constants. Defining 

w {a,x,y) = y a w [a,x,y) (149) 
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then gT: (x, y) and g™ (x, y) will satisfy the same ordinary differential equation 
if z (a, x, y) and w (a, x, y) satisfy the same pair of partial differential equations. 
With the notation 

dw „ dw „ d 2 w „ d 2 w - d 2 u> 
ox ay ctar oxoy oy z 

then 

w (a,x,y) = y~ a w (a,x,y) (150) 
P = jT Q p (151) 

q = y- a (q-~w) (152) 
r = y- a r (153) 
S = jT a f $ - ^ (154) 

t = y-*(i+^±^ w -2^) (155) 

and substituting these relations in equations (|145[) and (|146[) gives 

yt- xs + (1 - a) q + w = 0. (156) 

yu> = [2a — (2a + 1) x]p + a; (1 — x) f — y 2 i — a 2 w + (2a — 1) yq + 2xys. (157) 
Eliminating yw from (1157[) using equation (|156p gives 

x (1 — x) f + xys + [2a — (2a + 1) x]p + ay? — a 2 w = 0. (158) 

Equations (|156[) and (|158[) obtained from the Kampe de Feriet function are iden- 
tical to equations ([99]) and (|100|) obtained from the Horn function, so g™ (x, y) 
also satisfies the differential equation (| 1 33[) . 

5 Fourth-order differential equations from an in- 
tegral representation 

The integral representation (jTTJ) for the Fourier coefficient can be written in the 
form: 

(rR,\u) = (159) 



where 



. exp (iX-Jbj — cos ib) , , 
(A, w) = / Vy i — '- cos (rmb) dxb. (160) 
Vcj — cos ib 
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The function y m (A, u>) satisfies the partial differential equation 

d 2 y m {\uo) A 



dujdX 2 
and this has the elementary separated solution 



Vm (A, w) 



Y m (s, A, u>) = C m (s) exp 



A 2x 1 
±i [ su> H 

4s 



(161) 



(162) 



where s is the separation constant. The solution y m (A, w) can be constructed 
as a superposition of the allowable (i.e. finite at infinity) elementary solutions 
given by (|f 62|) . This gives y m (A,w) in the form: 



y m (X,ui) = / C m (s)exp 



A 

±i [ suj H 

4s 



(163) 



where the ± sign is chosen so that the integral converges. As u> is real and 
positive, this depends only on the imaginary part of A = a + ia, the appropriate 
sign being the same as that of cr, which will be assumed positive here. Setting 
A = in (fT60l) and (fT63]) gives 



(164) 



C m (s) exp (isuj) ds = V2Q m -i/2(v) 
and C m (s) can be determined from the integral [7, eqn 6.62f f] 

oo 

J ex-p(-s5)J„(s)s f *- 1 ds = 
o 

r {v + n) ^(v + nv + ii + i 



v + 1; =■ 

S 2 



(165) 



Setting S — —iu>, v — m and /i = 1/2 in equation (|165|) gives 

exp (iljs) J m (s) s -1 / 2 ds = 
r (to +1/2) /to +1/2 to + 3/2 



(_^m+l/2 2mbJ m+l/2 T ( m + 1) 



m + 1: 



(166) 



and the expression for Q m _i/2 (w) in terms of the Gauss hypergeometric func- 
tion is [16, eqn 8.1.3] 



i Fk r (to + 1/2) 

g m -i/ 2 h - y 2 ^+V2f (to + 1) 2 ^ 



771+1/2 TO + 3/2 1 

o ' o ' m + 1; ~~2 

2 2 w z 



(167) 
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From equations fTM]) . (fTTO)) and p^7]) it follows that 

c m ( s ) = V^H) m+1/2 J m (s) s - 1/2 

and hence 



(168) 



y rn (X,uj) = V^M)" 1 ^ 2 



exp 



4s 



J m (s) s-Wds. 



(169) 



For the special case of evanescent waves such that A = ia with a > then with 
a suitable transformation in the complex plane, this equation can be expressed 
in the form 



y m (ia,uj) = y/n exp 



4s 



I m (s) s-^ds. 



(170) 



The details of this transformation are given in Appendix B. Equation (|170p is 
straightforward to evaluate numerically as the integrand is not oscillatory and 
decays exponentially to zero as s — > oo provided u> > 1, which from equation 
(fT5|) is always the case. This follows immediately from the leading term in the 
asymptotic approximation as s — > oo of I m (s), which is [7, eqn 8.451 5]: 



I m (s) 



exp (s) 



(171) 



5.1 Fourth-order differential equation in terms of u) 

The integral representation (|169[) allows the ordinary differential equations in 
terms of w or A satisfied by y m (X,u>) to be constructed in a straightforward 
manner. It is convenient to define a new variable x such that 



A 



X 



and also a new dependent variable y m (x,^) such that 

Vm{X,U) = Vm (A,W) . 

Then y m (x,u) is given by 



(172) 



(173) 



J m (s)f(s,u,x)ds 



(174) 



where x is to be regarded as a constant embedded parameter in the ODE sat- 
isfied by y m (x,u), and where /(s,cj,%) is given by 



(175) 
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The various derivatives of y (a;, \) are then given by 

oo 

(if I ' s n J m (s)f(s,cj, X )ds. (176) 



d n ii 



duo 7 



i d ( s dJ 2 ^sy ={ ^__ 1]Jm{s) (177) 



The Bessel function J m (s) satisfies the differential equation [7, eqn 8.401] 

\ d ( dJ m (s) \ / 771* 
s V ds / \ s 2 

and therefore: 

oo 

' m 2 \ . . . fid/ dJ m (s) 



-l)J m (s)f(s,u, X )ds= I -—l 8 -^)f(8,cj,x)d8. (178) 







Differentiating (I178|) twice with respect to w and utilizing equations (I174|) - ()176|) 
gives: 

oo 

(m 2 - s 2 ) J m (s) f(s, u), x)ds = J s-jj- ^ S ~^~^ J Z( s ' w ' x)ds 

o 

oo 

+ m y m = / s— s — f(s,u),x)ds. (179) 



cZw 2 a '" J ds \ ds 

o 

Integrating this twice by parts yields 

oo 

^ + m 2 y m = J J m ( S )± (s^ [ S /(.s_v X )] j ,/,, (ISO) 


Since 

Ts { s Ts [s/(s ' w ' x)] ) = (l + 2 * ws " wV + 2wx ~ x2s_2 ) /(s ' w ' x) (181) 

then 

oo 

~^JT + ?™ 2 y m = J J m{s) Q + 2iws - to 2 s 2 + 2ujx - X 2 s~ 2 ^j f{s, w, x)ds. 

(182) 

For the Poisson case with x = 0, employing (|176p in (|182|) gives Legendre's 
equation (pD|) of degree m — 1/2, as must be the case for consistency with 
equation (fl"4"|) . For the Helmholtz case, (|182[) must be differentiated twice with 
respect to ta before employing (|1T6[) . This yields the fourth-order linear ODE 

(1 _ ^ 2) <H^ _ + (m2 _ ^d^_ _ ^ = Q (183) 

which becomes equation (1551) on substituting x = A 2 /4. 
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5.2 Equivalence of the differential equations 

Equation (|183p can be converted to an equation in terms of x = k 2 by making 
the substitutions: 

X = f (184) 



2-x 



x 



in 



(185) 



1-^ = 11^1) (186) 

x z 

25 

- 2 X uj - — = a 2 - a- 6+ {x-2)y (187) 

dy x 2 dy 
dm 2 dx 
d 2 y x A d 2 y x 3 dy 
duj 2 4 dx 2 2 dx 
d 3 y a; 6 d 3 y 3a; 5 cPy 3a; 4 
duj 3 ~~ 8 dx 3 4 dx 2 4 efe 
d 4 y _ a; 8 d 4 y 3x 7 d 3 y 9x 6 d 2 y 3a; 5 dy 
d^ = I6dx 1+ ~dx^^ ~dx 2+ ~dx' 
On collecting terms and simplifying this yields 

(l-a;)a; 4 + (6-9x)a; 3 + 

d 2 v 

\a - a 2 + 6 - 18a; + (2 - x) y] x 2 —^- + 
1 aar 

-2[a(a~l) + a;(3 + 2/)]a;^+2/ 2 y = 0. (192) 

dx 

Inspection of equations (|133p and (|192[) . obtained by totally different methods, 
shows that they are identical. 



(188) 
(189) 
(190) 
(191) 



6 Comments and conclusions 

The Fourier coefficients for the Helmholtz Green function have been split into 
their half advanced+half retarded and half advanced— half retarded compo- 
nents, and these components have been given in closed form in terms of two- 
dimensional hypergeometric functions. These solutions generalize the well- 
known solutions of Poisson's equation for ring sources, and reduce to them 
in the static limit when the wave number f3 = 0. The two-dimensional hy- 
pergeometric functions can be considered as double series, with the order of 
summation arbitrary. The two summation choices give different series of special 
functions for each of the Fourier components, and all of these series have been 
numerically verified, as have the closed form solutions themselves. One series is 
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given in terms of Hankcl functions, and only a few terms are need far from the 
ring source for accurate results. A second series in terms of associated Legen- 
dre functions only requires a few terms in the neigborhood of the ring to give 
accurate results. 

The systems of partial differential equations associated with each of the two 
generalized hypergeometric functions have been used to derive a fourth-order 
ordinary differential equation in terms of x = k 2 for the Fourier coefficients. A 
completely different approach involving integral representations of the Fourier 
coefficients has been presented in tandem, which derives many of the same 
results, as well as some new ones. Both approaches give exactly the same fourth- 
order differential equation for the general Fourier coefficient, despite the algebra 
being rather intricate in both cases. Another fourth order ordinary differential 
equation in terms of the wave number parameter A can also be derived by the 
methods presented here. 

A Series from an integral representation 

The Fourier coefficient (/?, r, R, z — Z) = G r ^ (rR, A, uo) given by equations 
(|27|) and ([28]) can be expressed in the form 



G£ {rR, A, w) = j J2 [ {to - cos t/.) (, ^ 1)/2 cos (mip) (193) 

77 rR ™=° o n ' 

where u> and A are defined by equations (|15|) and (|43| respectively. Evaluation 
of (|193p requires the integral 

im,n{ L0 )= / cos (mip) (uj — cos ip)^" 1 ^ 2 dtp (194) 





which can be evaluated for m G No using the integral representation [7, eqn 
8.711 2]: 

p - m(g) = r lr|^t) 1) / cos (^)(e + v / ^ cos0 ) l/ ^ ( 195 ) 



which is equivalent to 

(-i)-r(i/+m+i)(e 2 -ir /2 7 / f y 

P " l ^= WtTT ycosM^^-cos^ # . 

(196) 

The substitutions 

u = (197) 
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and 



then give 



n-1 



(-l)"Vr((n+l)/2) 

r (m + (n + 1) /2) 



(n-l)/4 



(n-l)/2 



(198) 



(199) 



and (|193[) becomes 



(-If 



E 

n=0 



(*A(, 



c 2 -l) 1/4 



(w 2 - 1) 1/4 V2^i? 
r((n+l)/2) 



pn 
1 1 , 



i! r(m+(n + l)/2) ( "- 1)/2 Va/^ 2- T 



(200) 



Splitting the series (|200[) into even and odd terms gives after some reduction 

(-irvsF „ 



(w 2 - 1) 1/4 7277?' 



E 

p=0 



p." 



pW(p + m + l/2) p - 1/2 \V^l 



(201) 



A™M,A lW ) = H^x 



2V2rR 



E 



A\/i?-T 



r(s + 3/2)r(m + s + i) s VV7 2- T 



where the factorial formulas 



(2p) 



! — 



2 2 *T(p + l/2)p! 



(2s + 1)! = 



2 2s+1 r(s + 3/2)s! 



(202) 



(203) 



(204) 



have been used to obtain (|20ip and (|202|) . The index in (|202p runs from s = m 
rather than from s = as the associated Legendre polynomial P" 1 (£) is zero 
for s < m. The substitution s = p + m in (|202p gives the alternative form 



'in 



A X 2 yAj T ^l 



E 

p=0 



-A 2 ^^ 2- 1 



2V2rR 
1 



r(p + m + 3/2)r(p + 2m + l) p+m V V^~T J ' 



(205) 
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The indices in equations (|201[) and (|205|) can be switched to negative values 
using the relations [16, eqns 8.2.5, 8.2.1]: 

1 (u — m + 1) 
PZ-i (0 = (0 (207) 

which gives 

A" 1 {rR, A, w) = ^ttt^ x 

+ (c; 2 -l) 1/4 V2^R 



4 J n \v(n- m , + \ii\ p -p-vz ( .ra-r r) (208) 



p=0 



A™ (rfl,A,w) = 



pir(p-m+l/2) -p-^Va/^ - T 

y^TTlX I X 2 y/u) 2 — 1 



2V277? 



-AV^^iV i 



4 J r(p + m + 3/2)r( ? ,+ i) P -™-™- 1 ( v vZ^Tj- (209) 

The kind of Legendre functions in equations l|208p and (|209p can be switched 
using the Whipple relation [16, eqn 8. 2. 7], [17]: 



a _ 1/a / ^ \ _ (^ 2 -l) 1/4 expH M 7r) 

-"- 1/a I V^i ) ~ (V2) 1/2 r (a + „ + 1) Qq { ) ( } 



which gives after some reduction 



A™(ri?,A,^) = ^Lx 



f WZEl Y gizVjM (211) 

/-^ 4 nir — m 4- 1 r fry) 4- n 4- 1 1 0\ V ' 



p=0 



p\T (p - m + 1/2) T (m + p + 1/2) 
f-l) m 

A m (ri?,A,w) = v Z__ x 



^ AW^l Qm-1/2 M (21) 

2—1 A I nir (n 4- -m -I- .1 /0\ r (r> 4- 9<m 4- 1 V 1 ' 



p=0 
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B Transformation of an integral in the complex 
plane 

The integral for y m (ia, to) in equation (|169p for A = ia can be considered to be 
the contribution along the real axis of the contour integral 

J m (s)s- 1/2 ds (213) 
c 

where C is the closed contour shown in the figure below, in the limits as e — > 
and R — > oo. 



m+1/2 



exp 



a 
is 



Figure 1: Contour C in the complex s— plane where s = x + iy. The smaller 
quarter circle is of radius e — > and the larger quarter circle is of radius R — > oo. 
The integrand of the contour integral has an isolated essential singularity at the 
origin of the s— plane and is analytic on and within the contour C . 

The corresponding contribution to the contour integral along the imaginary 
axis is given by 



[/ m ( U ,a) = -^Hr +1/2 



exp 



uiy - 



Ml 



J m {iy)y- 1/2 W 1/2 rfy 



(214) 

and this can be stated in terms of the modified Bessel function of the first kind 
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using the identity [7, eqn 8.406 3] 

J m (iy) = (i) m i m (y) 

which gives immediately 



°° 2 

U m (cr,ui) = J cxp -^ujy+^j 



(215) 



I m (y)y- 1 / 2 dy. (216) 



The integrand of (|213[) is analytic everywhere within the contour C and therefore 
from Cauchy's theorem, the integral (|213[) is zero. Therefore if the contributions 
to the contour integral along the two quarter circles vanish in the limits as e — > 
and R — > oo, then 

Vm ( icr , <*>) + U m (cr, u) = (217) 

and equation p70[) has been proven. Along the smaller quarter circle we set 
s = e exp (iff) and the contribution to the contour integral becomes 



tt/2 



V m (a,cj) = -V^(-i) m+1/2 



exp 



i ( we exp (iff) — — exp (—id 



J m (eeyLv(i6))ie 1/2 d9 



(218) 



and this clearly vanishes as e — > 0. Along the larger quarter circle we set 
s = Rexp (iff) and the contribution to the integral is given by 



-it/2 

W m (a,w) = s/tt (-i) m+1/2 J exp i y 



i ( LuRexjp (iff) — exp (—it 



J m (R exp (iff)) iR 1/2 exp (i0/2) d0. 



(219) 



The leading term in the asymptotic approximation of J m (s) is given by [7, eqn 
8.451 1] 



J m (s) 



2 / wm 7r s 

cos s 

its V 2 4> 



(220) 



and therefore 



tt/2 



W m (a,u) = ^{-iT +1/2 
/•\l/2 r 

(0 



exp 



i ( LuRexp (iff) — exp (— it 



o 



exp ii? exp (id 



mm iit 
~2 T 



exp ( — ii? exp (i0) 



i7rm in 



d6 
(221) 
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which can be expressed as 



exp 



+ l)Rexp(i6 



7T/2 



exp \ -i—cxp(-tb 



~~2 T 



exp 



i ( (u - 1) R exp (z0) + — ^— + 
(222) 



Inspection of (|222p shows that the integrand vanishes exponentially in the limit 
as R — > oo provided u) > 1. This condition always holds and is also the condition 
for the integral in (|170[) to converge. 



C Numerical results 

The series solutions for the Fourier coefficients given by equation (|54p and equa- 
tions (jT?l) - (|T8")) were evaluated using Mathematica and the numerical perfor- 
mance was explored for various geometric parameters and wave numbers. For 
comparison, the two integrals (fTTj) and (|2"2"|) for the Fourier coefficients were also 
evaluated numerically for the same parameters. All four methods give iden- 
tical results at locations which are neither too far away nor too close to the 
ring source. The numerical integration (fTTj) performs very well at all distances 
from the ring, whereas the numerical integration (|22p fails when either very 
close to the ring or too far away. No cases were identified where equation (|2"2"j) 
was superior. The Hankel function series requires fewer and fewer terms 
for convergence as the distance from the loop increases, and conversely perfor- 
mance decreases as the ring is approached. The associated Legendre function 
series ([TT]) and (fT5)) have precisely the opposite performance, with great accu- 
racy close to the ring and failure at large distances from the ring. The two series 
(|4"T|) and (|4"5|) are well suited to calculations close to the ring as the associated 
Legendre functions themselves each contain the ring singularity as to — > 1. By 
contrast, the Hankel functions (|54j) are not singular at the ring and hence an 
increasing number of terms are required to model the singularity as the ring is 
approached. In all cases, there is are always at least one numerical integration 
and one series solution which can be used to cross check each other. 

Sample numerical results are given in Table 2 for moderate distances from 
the ring source, and shows the number of terms required by each series to 
match the numerical integrations exactly. Table 3 shows the performance of the 
Hankel series with increasing distance from the ring. The number of Hankel 
terms decreases to very few at large distances from the ring. Table 4 shows 
the performance of the two associated Legendre series ([4"T]) and (|4"5j) as the ring 
is approached. The real part of G^ ((3, r,R,z — Z) diverges logarithmically as 
z — > Z, whereas the imaginary part tends to a finite limit. It can be seen 
immediately from Table 4 that only 8 terms in each Legendre series is sufficient 
to calculate ((3, r,R,z — Z) for the range {0 < z ^ 1}. 



31 



m 


(3 


r 


z 


N t 


N 2 


G%((3,r,R,z-Z) 





2 


1/2 


1/2 


94 


10 


-0.4332208795 + 0.6063507453 i 





2 


1/2 


3/2 


24 


12 


-0.4324244083- 0.2593676946 i 





2 


1/2 


5 


9 


24 


-0.1320141290 - 0.1413052175 i 





2 


1/2 


10 


8 


37 


0.02892833221 + 0.09482283693 i 





2 


1/2 


20 


6 


63 


-0.03552779935 + 0.03502697525 i 


3 


5 


3/2 


1/2 


224 


24 


-0.2152817201 - 0.2085849956 i 


3 


5 


3/2 


1 


97 


25 


0.1382226177- 0.2010980843 i 


3 


5 


3/2 


5 


17 


47 


-0.009794158906- 0.000546039281 i 


3 


5 


3/2 


10 


13 


80 


-0.0004846328044+0.0006340532520 i 


3 


5 


3/2 


20 


10 


147 


0.00000356468968 + 0.00005347913049 i 



Table 2: Numerical solution for the Fourier coefficient G^ (/3, r, R, z — Z). R — 
1, Z = and the other parameters are given in the table. N\ is the number 
of terms in the Hankel function series (|54[) needed to give the accuracy given. 
N2 is the number of terms in the associated Legendre series (|47p and (|4"5|) to 
provide the accuracy given. Of the two associated Legendre series, (|47l) requires 
more terms than (14"51) for the stated accuracy. 



z 


Ni 


G^ (/J, r 


R,z-Z) 


1/2 


227 


0.0785417676- 


- 0.2281496125 i 


1 


99 


0.1318397799- 


h 0.0959755332 i 


5 


20 


0.04717552085- 


- 0.09819984770 i 


50 


7 


-0.001762722093 


- 0.000313668126 i 


100 


6 


-0.0000246835014 


+ 0.0004487208692 i 


200 


5 


-4.36384289 x 10~ 6 


- 0.00011237773355 i 


1000 


4 


-1.884281670 x 10~ 6 


- 4.086433833 x 10" 6 i 


5000 


3 


-1.446923432 x 10" 7 


+ 1.070704963 x 10~ 7 i 


10000 


2 


4.307310056 x 10~ 8 - 


h 1.302718185 x 10~ 8 i 


10 7 


1 


-2.303180610 x 10" 14 


+ 3.865922798 x 10~ 14 i 



Table 3: Series solution for the Fourier coefficient G^ (j3, r, R, z — Z). m = 1, 
(3 — 6, R=l, Z = Q and r = 3/2. z is given in the table. N± is the number of 
terms in the Hankel function series (|54p needed for the accuracy given. 
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z 






Cr H r 


, il, 2. — Z/ J 


1 




8 


0.1874175169- 


f 0.1222388714 i 


10- 


1 


8 


0.8955546890- 


f 0.1360159497 z 


10" 


2 


8 


1.628566013- 


f 0.136158894 i 


10" 


3 


7 


2.361506874- 


f 0.136160324 i 


io- 


4 


7 


3.094442571 - 


f 0.136160339 i 


io- 


5 


7 


3.827378171- 


f 0.136160339 i 


io- 


6 


7 


4.560313770- 


f 0.136160339 i 


io- 


7 


7 


5.293249369- 


f 0.136160339 i 


10" 


8 


7 


6.026184968- 


f 0.136160339 i 


10" 


9 


7 


6.759120567- 


f 0.136160339 i 



Table 4: Series solution for the Fourier coefficient (/?, r,R, z — Z). m = 2. 
(3 = 1, R = l, Z = and r — 1. z is given in the table, and as z — > 0, the ring 
is approached. A^2 is the number of terms in the associated Legendre function 
series (|47|) and (|48[) needed for the accuracy given. 
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